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ABSTRACT 

We present high-contrast images of HR 4796 A taken with Subaru/HiCIAO in 77-band, resolving the debris disk 
in scattered light. The application of specialized angular differential imaging methods (ADI) allows us to trace 
the inner edge of the disk with high precision, and reveals a pair of “streamers” extending radially outwards 
from the ansae. Using a simple disk model with a power-law surface brightness profile, we demonstrate that the 
observed streamers can be understood as part of the smoothly tapered outer boundary of the debris disk, which 
is most visible at the ansae. Our observations are consistent with the expected result of a narrow planetesimal 
ring being ground up in a collisional cascade, yielding dust with a wide range of grain sizes. Radiation forces 
leave large grains in the ring and push smaller grains onto elliptical, or even hyperbolic trajectories. We measure 
and characterize the disk’s surface brightness profile, and confirm the previously suspected offset of the disk’s 
center from the star’s position along the ring’s major axis. Furthermore, we present first evidence for an offset 
along the minor axis. Such offsets are commonly viewed as signposts for the presence of unseen planets within 
a disk’s cavity. Our images also offer new constraints on the presence of companions down to the planetary 
mass regime (~9 Af Jup at 075, ~3 A/j up at 1"). 

Subject headings: circumstellar matter — planetary systems — techniques: high angular resolution - stars: 
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1. INTRODUCTION 

Vega-type debris disks were first identified by infrared (IR) 
excesses around nearby main sequence stars { Aumann et al. 
1 984). The dust content of these second-generation disk sys- 
tems is believed to be continuously replenished via collisional 
breakup of remnant planetesimals (cf. Wyatt 2008). 

Since the imaging of the 5 Pic system (Smith & Terrile 
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1984), nearly two dozen nearby debris disks have been 
spatially resolved. The morphological appearance of re- 
solved debris disks is predicted to be influenced by in- 
teractions between dust in the disk and nearby planets 
(Ozemoyetal. 2000; Kuchner & Holman 2003), the local 
interstellar medium (Maness et al. 2009), recent stellar fly- 
bys and binary companions (Wyatt et al. 1999), mutual grain 
collisions (Thebault& Wu 2008), and interaction of dust 
with residual gas (Klahr & Lin 2005). Many resolved 
systems exhibit all morphological structures predicted by 
these mechanisms (Schneider etal. 1999; Kalas et al. 2005; 
Golimowski et al. 2006; Kalas et al. 2008). The observable 
morphology of resolved, optically thin debris disks is also 
wavelength dependent, as different bandpasses sample differ- 
ent grain size distributions (Wyatt 2006). 

HR 4796 A is a young (~8-10Myr; Stauffer et al. 1995), 
nearby (72.8 ±1.7 pc; van Leeuwen 2007), AOV-type star 
first identified as a debris disk system from an IR excess 
observed with IRAS (Jura 1991). It has a co-moving M-type 
stellar companion at a separation of l".l (Jura etal. 1993). 
The HR 4796 A debris disk has been spatially resolved 
at numerous optical, infrared, and sub-millimeter wave- 
lengths (e.g., Koemeretal. 1998; Jayawardhana et al. 1998; 
Schneider et al. 1999; Sheret et al. 2004; Hinkley et al. 2009; 
Schneider et al. 2009), appearing as a narrow, highly inclined 
belt with a radius of 1705. Models suggest a dual dust popu- 
lation including gravitationally confined grains at ~9 - 1 0 AU 
in addition to those imaged at ~70 AU (Augereau et al. 
1999), possibly involving radiationally evolved organic 
materials (Debesetal, 2008). Wahhajetal. (2005) present 
multi-wavelength observations showing that the outer disk 
cannot be modeled as a single dust component; these authors 
adopt a two-component model with two different grain sizes. 
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Fig. 1 . — High-contrast images of the HR 4796 A debris disk. White plus signs mark the star's location. The saturated central area is masked in grey, (a) 
Subaru HiClAO 77 -band data after reduction with simple ADI (subtraction of a median background) and median smoothing on the scale of 5 pixels (= 48 mas 
ss 1 A/D). Orange dotted lines mark the mean of the spider diffraction pattern. The stretch is linear between ±3- 10“ B times the stellar peak flux, (b) The same 
after ADI reduction with the LOCI algorithm, using conservative settings (optimization region area of 10.000 PSF footprints). Note the “streamers” extending 
radially out from the ansae (yellow arrows). Linear stretch of ±1.5 • 10” 5 . (c) The same with aggressive LOCI settings (300 PSF footprints) for point-source 
detection. Linear scale of ±6 • 10 ” 6 . (d) HST/STIS 0.2-1.0/j.m image front Schneider et al. (2009, Figure 8) reprinted for comparison. These authors point 
out a possible flux deficit in the southwestern quadrant of the disk (dashed green arrows), (e) Simple ADI image from (a) in logarithmic stretch [10~ 7 , 1 0 — 3 ] , 
showcasing the streamers. Red shading indicates spider artifacts. 


Here we present new ground-based near-infrared imaging 
data at high spatial resolution, revealing the tapered outer 
regions of the debris disk. 

2. OBSERVATIONS AND DATA REDUCTION 

Our observations of HR 4796 A were obtained with the 
Subaru Telescope on May 24, 2011, within the SEEDS sur- 
vey (Strategic Exploration of Exoplanets and Disks with 
Subaru'HiCIAO, Tamura 2009). The HiCIAO instrument 
(Hodapp et al. 2008) with a 20" x 20" field of view and a plate 
scale of 9.50±0.02mas/pixel was used. The image rotator 
operated in pupil-tracking mode to enable angular differential 
imaging (ADI, Marois et al. 2006). A sequence of 260 images 
was taken in II -band with an exposure time of 10 s, for a total 
integration time of 43.3 min and a total field rotation of 23°. 
Weather conditions were good (seeing 0 / /5—0 , /8 in U-band), 
and the A0188 adaptive optics system (Minowaetal. 2010) 
provided a FWHM of 6.5 pixels = 62 mas. 

The images were corrected for Hatfield and field distortion 
(Suzuki et al. 2010). Stellar position was estimated in each 
frame by two different methods: (1) Fitting a Moffat profile 
to the PSF halo, and (2) triangulating between symmetrical 
pairs of static speckles. Both were consistent with an empir- 
ical drift model (r, = f*o + vo'i + 1/2 do* 2 ) plus measure- 
ment noise. The difference between the two sets of centroids 
showed no systematic behavior, and was consistent with in- 
coherent combination of the two measurement noise sources. 
We therefore used the drift model for image registration, for 
an estimated centering accuracy of ~0.3 pixels = 2.9 mas in 
the co-added image. 

ADI combined with the LOCI algorithm (Locally Op- 
timized Combination of Images, Lafreniere et al, 2007) 
is currently the most successful ground-based imag- 
ing technique for the detection of planets (Marois et al. 
2010; Lagrange et al. 2010; Currie et al. 2011) and substel- 
lar companions (Thalmann et al. 2009; Biller et al. 2010; 
Jansonetal. 2011). Furthermore, it has proven useful in 
revealing hiah-contrast circumstellar disks (Thalmann et al. 
2010: Buenzli et al. 2010). 

We applied three implementations of the ADI technique to 
our data: (I) “Simple ADI"’, consisting of subtracting a me- 


dian background from the entire dataset before derotating and 
co-adding. This method causes the least amount of flux loss 
and is therefore useful for estimating the surface brightness 
profile of the disk. (II) “Aggressive LOCI”, using frame se- 
lection criteria of 0.5 FWHM (minimum differential field ro- 
tation between images to be used for mutual subtraction, to 
limit self-subtraction of physical sources) and an optimiza- 
tion region with an area of 300 PSF footprints. Although most 
of the disk signal is lost and negative over-subtraction effects 
appear, this method achieves excellent speckle suppression 
and provides the best constraints on point sources, such as 
planets. (Ill) “Conservative LOCI”, a compromise between 
the previous two methods, first described in Thalmann et al. 
(2010). This method preserves more disk flux than aggressive 
LOCI while achieving significantly better speckle subtraction 
than simple ADI. The resulting image proves ideal for deriv- 
ing geometric parameters of the inner disk edge. Due to the 
dataset’s limited field rotation budget, we do not use an in- 
creased frame selection criterion to conserve more disk flux, 
as Thalmann et al. (2010) do. Instead, we enlarge the opti- 
mization area to 10,000 PSF footprints, lowering the impact 
of the disk signal on the optimization process, as first demon- 
strated in Buenzli et al. (2010). Unless otherwise noted, we 
use the default numerical and geometric LOCI settings as de- 
fined in Table 1 and Fig. 1 of Lafreniere et al. (2007). 

3. RESULTS 

3.1. Imaging the debris disk 

Figure 1 shows the results of applying the three ADI re- 
duction methods to our data: Simple ADI in (a), conservative 
LOCI in (b), and aggressive LOCI in (c). In panel (d), we 
reprint the HST/STIS image by Schneider et al. (2009, Figure 
8) for comparison. Starlight scattered from the disk is visible 
as a strongly projected ellipse in all images, down to stelloeen- 
tric separations of 07 5 in simple ADI and 074 in LOCI. Dif- 
ferential imaging clearly resolves the inner edge of the disk. 

While most of the flux concentrates around this edge, each 
ansa appears radially smeared out, extending a “streamer” 
outwards along the major axis. Similar streamers have been 
observed in deep ADI or roll subtraction images of other de- 
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Extended Outer Regions of the Debris Ring Around HR 4796 A 


TABLE l 

Numerical results on HR 4796 A 


Maximum merit ellipse size 

ang. sep. 

projected (AU)* 1 ) 

Semimajor axis uW 

17087 ± 07023 

79.2 ± 1.7 

Semiminor axis b 

07248 ± 07 007 

18.0 ±0.5 

Maximum merit ellipse orientation 

degrees 


Position angle U 

26.4 ± 0.5 


Inclination i — arccos (b/a) 

76.7 ±0.5 


Maximum merit ellipse center 

ang. sep. (mas) 

projected (AU)* 1 ' 

x (along RA) 

-6.6 ± 4.0 (3) 

-0.48 ±0.21® 

y (along Dec) 

-22.1 ±4.8< 3 > 

-1.61 ± 0.28* 3 ' 

u (along major axis)* 4 ' 

— 16.9 ± 5.1 

-1.23 ±0.31 

v (along minor axis)* 4 ' 

15.8 ±3.6 

1.15 ±0.16 

Upper limits on planets 

5 a contrast 

Mass (M Jw )* 5 ' 

At CK'3 ==18 AU in proj. 

T 

o 

oo 

~17 

At O'/S = 36 AU in proj. 

■j 

f 

o 

~9 

At 1"= 73 AU in proj. 

5.7- I0“ 6 

~2.8 

At 2"— 146 AU in proj. 

8.6 ■ 10- 7 

~1.4 


NOTES. (1) Conversion from arcsec to projected AU assuming a distance of 
72.8 pc; the uncertainty of 1 .7 pc is not included in the errors. (2) Based on the 
model disk described in Section 4. 1 , the photometric peak-to-peak radius a is 
~96% of the maximum merit semimajor axis; this would yield a value of a = 
l'/044 ± 0.023 consistent with Schneider et al. (2009). (3) Errors in ( x , y) 
are correlated; use the uncorrelated errors in (u, v) instead. (4) The axes 
of the (u, v) coordinate system are aligned with the disk’s major and minor 
axes, thus they are rotated counterclockwise from the RA/Dec coordinate 
system (x, y) by Q = 26.4°. (5) Conversion from //-band contrast via the 
COND evolutionary models (Allard et al. 2001 ; Baraffe et al. 2003). assuming 
an age of 8 Myr. 


bris disks, such as HD 61005 (“The Moth”, Buenzlietal. 
2010) and Fomalhaut (Kalas et al. 2010). They are thought 
to represent part of the thin, extended, tapered outer reaches 
of the debris disk, which are most detectable at the ansae. 
Oversubtraction makes them appear narrower than they are. 
In both simple ADI and LOCI, the flux at a given separation 
is unaffected by any flux at smaller separations (Marois et al. 
2006; Lafreniere et al. 2007), thus the streamers cannot be ar- 
tifacts caused by the bright ring further inwards. 

Schneider et al. (2009) note a marginal flux deficit in the 
southwestern quadrant of the debris disk. Although there are 
some indications of substructure in the same location in our 
dataset (dashed green arrows), the S/N ratio is insufficient to 
confirm the existence of a physical disturbance in the disk. 

3.2. Ring geometry from maximum regional merit fitting 

The crisp representation of the disk’s inner edge in the con- 
servative LOCI image allows accurate measurement of its ge- 
ometrical properties. For that purpose, we have devised the 
maximum regional merit technique, which is loosely based on 
the method of Buenzli et al. (2010). 

First, we apply a median filter on the spatial scale of 5 pix- 
els («1 A/D), which reduces the pixel-to-pixel noise while 
keeping larger edges and structures intact. We then convert 
the image into a signal-to-noise map (S/N; Figure 2a) by cal- 
culating the standard deviation in concentric annuli around the 
star and dividing the pixel values in each annulus by this noise 
level. The disk itself is masked during noise calculation. The 
S/N map emphasizes the well-resolved outer parts of the disk. 

We then generate a large number of ellipses with random- 
ized major and minor semi-axes (a.b), position angles (fi), 
and center positions relative to the star (x, y). We assign a 
merit value r; to each set of parameters, defined as the mean 
value of the pixels in the S/N map that lie within an ellipse- 
shaped sampling area described by the numerical parameters. 
The sampling area for a given parameter set (a, b, SI, x, y) is 





-20 -10 0 10 
East/West offset (mas) 


Fig. 2. — Derivation of geometric parameters, (a) Signal-to-noise map of 
the conservative LOCI image of HR 4796 A, calculated in concentric annuli 
around the star and excluding the bright ring from noise estimation. The 
stretch is [—10 cr, 10 or), (b) The same, but with red shading marking the 
sampling ellipse pixels for the best-estimate ellipse, (c) Histogram of merits 
rj measured in empty sky, showing a quasi-Gaussian distribution, (d) Scatter 
plot of the center offset parameters (x, y) for the solution with merits rj within 
1 <Jr) (black), 2 or) (red), and 3 u n (orange) of the maximum merit value T^ax, 
drawn from Monte Carlo simulations of 300,000 parameter sets. A white plus 
sign marks location of the best estimate (5, y): a black plus sign marks that 
of the star. To guide the eye, white dashed ellipses with semiaxes of 1-4 
spatial standard deviations of the 1 <j fl family along the major and minor axes 
are provided, as well as dotted black circles marking the 1-4 a uncertainty 
regions of the star’s position. 


defined as the area between two nested ellipses that share the 
same orientation and center (O, x , y) but have semi-axes of 
(a + 6, b + e) and (a — S,b — e), respectively. We choose 
£ = 0.7 px, which results in a single-pixel thickness along the 
minor axis, and S — 4px, for a 7-pixel thickness along the 
major axis. The value for <5 corresponds to the observed width 
of the high-flux areas in the ansae of the disk; this is to ensure 
that a well-fitted sampling ellipse will be delimited by strong 
gradients in the S/N image on all sides, and thus sensitive to 
displacement and scaling (Fig. 2b). 

We find that this method yields more accurate results for 
our data than the conventional method of fitting of an el- 
lipse to flux maxima in the ring (as applied to HST data by 
Schneider et al. 2009). since our data has finer spatial sam- 
pling but a higher photometric noise level than the HST data. 
The disk profile is flat around the flux maxima and thus vul- 
nerable to noise, while the tracking of sharp gradients is ro- 
bust and profits from the fine spatial sampling. 

We find a maximum merit of r) mm = 3.46. In order to 
estimate the noise inherent to this merit figure, we measure 
n for 10 4 sampling ellipses placed in the empty sky region 
90° away from the position angle of the disk. W’e find a 
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Fig. 3 - Modeling of the HR 4796 A debris disk, (a) Simulated scattered- 
light image of the simple power-law model disk described in the text, dis- 
played in linear stretch, (b) The same in logarithmic stretch, (c) PSF- 
subtracted HST/ST1S image from Schneider et al. (2009) in logarithmic 
stretch. Note that a faint halo as predicted by the model is visible beyond 
the bright ring, (d) Result of conservative LOCI reduction of our Subaru 
HiCIAO data after injection of the model disk at a 90° position angle off- 
set (linear scale). The model appears qualitatively indistinguishable from the 
real disk, demonstrating that the morphology of the disk signal in our ADI 
images is well understood. 

roughly Gaussian distribution (Fig. 2c) with a standard de- 
viation a n = 0.077 and a mean value p v = -0.0006 (« 0 
within the errors). We therefore define the family of “well- 
fitting” solutions as those satisfying rj > '/? max — a v , and the 
“best estimate” (a,b,Cl.x,y) as the mean value of the well- 
fitting solution family. 

Figure 2d plots the spread of the center offset values (x, y ) 
of solutions within 1 a Tj , 2a v , and 3 cr n of the maximum 
merit value ry*x- The contours of those solution families are 
concentric ellipses with axes aligned along the debris ring’s 
axes, demonstrating approximately Gaussian error behavior. 
To further test the validity of this maximum regional merit 
technique, we produced 10 model disks (as discussed in Sec- 
tion 4.1) with random offsets from the star, injected them into 
the empty sky region in our data, reduced them with conser- 
vative LOCI, and retrieved the ellipse parameters as described 


above. Deviations between measured and real offset coordi- 
nates were consistently below 1 a, with an RMS of 0.4 a. 

The best-estimate parameters and their errors, including 
both fitting errors and systematic errors from image registra- 
tion and plate scale, are presented in Table 1 . Most notably, 
we find that the ellipse center is offset along its major axis 
(towards the southwest) by u = 16.9 ± 5.1 mas, deviating 
from the star’s position at 3.3 a significance. Since the ma- 
jor axis is unaffected by foreshortening from projection along 
the line of sight, this translates directly into a physical dis- 
tance of 1.2 ± 0.3 AU, assuming an approximately circular 
cavity. Combined with Schneider et al. (2009)’s independent 
measurement of 19 ± 6 mas (3.2 a confidence), we can now 
confirm that HR 4796 A’s debris disk is offset from its star. 
Furthermore, our data suggest an offset along the minor axis 
of 1 5.8 ±3.6 mas in projection. If confirmed, this would cor- 
respond to a physical offset of ~5 AU. 

3.3. Constraints on planets 

Aggressive LOCI reduction of our image data does not 
detect any planet candidates in the system, but imposes 
upper limits on the masses of potential unseen planets. 
We use the photometric calibration method described e.g. 
in Thalmann et al. (2011), estimating and correcting the 
separation-dependent partial flux loss from LOCI processing 
as described in Lafreniere et al. (2007), and converting the H- 
band flux into planet mass using the COND evolutionary mod- 
els (Allard et al. 2001; Baraffe et al. 2003). Table 1 lists the 
resulting 5 o contrasts and planet masses for four separations. 
The noise is calculated in concentric annuli, and the spider 
artifacts are masked out for this purpose. 

4. MODELING 

4. 1 . Simple power-law model 

In order to understand the morphology of our disk images 
and gauge their usefulness for deriving disk properties, we 
generate a model disk by calculating synthetic scattered light 
images of inclined rings, with a power-law distribution of the 
scattering cross section R. We assume the inclination of the 
disk to be given by the axis ratio, and construct a Keplerian 
disk with a small but nonzero eccentricity to match our ob- 
served semi-major axis offset (1.2 AU, Table 1). The main 
free parameters of the disk model are the radial power-law 
slopes, 



where r is the radial distance to the star projected onto an 
ellipse with the desired geometry, «o is the semi-major axis, 
a m > 0 and ow < 0, and R(r) is related to the surface den- 
sity by (7 (r) = cr o • R(r) ■ (r /ao). We use a Henyey-Greenstein 
phase function and assume it is the same throughout the disk. 
To compute images, we use the single-scattering radiative 
transfer code GRaTeR (Augereau et al. 1 999) to compute syn- 
thetic images. This approach is valid for an optically thin disk. 

Given the scope of this model, we do not attempt to fit but 
merely approximate the real debris disk. Our working pa- 
rameters are a m = 35, indicating a very steep inner edge, 
Oout = —10, and an asymmetry factor g = 0.10, implying a 
slight preference for forward scattering (cf. g — 0.16 ± 0.06 
in Schneider et al. 2009; g = 0.03-0.06 in Debes et al. 2008). 
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Fig. 4. — Surface brightness profiles of the HR 4976 A debris disk (black 
curves, dotted 1 cr error bounds) extracted from the ADI image in Figure 3 
in strips 1 1 pixels wide oriented along the major axis. The blue curves mark 
the corresponding profiles in the model disk in the same image. The latter 
demonstrate that the straight slope of the injected model disk is preserved in 
simple ADI reduction. The red lines are linear fits to the plots in the interval 
If 1-1" 6 with respect to the disk center. The numbers indicate their slopes 
(i.e., the power-law exponents), including fitting errors. 


Figures 3a and 3b show the simulated scattered-light im- 
ages of the model disk at linear and logarithmic intensity 
scale, respectively. In Figure 3c, we show the PSF-subtracted 
image of the HST/STIS data from Schneider et al. (2009) in 
logarithmic scale for comparison, revealing an extended halo 
beyond the bright ring that matches the expected morphology. 
Since this dataset is reduced with PSF subtraction rather than 
ADI, it suffers from virtually no flux loss, but the background 
is dominated by residual speckle noise. 

To investigate the influence of ADI on disk flux, we inject 
the simulated scattered-light images of the model disk into 
our raw Subaru/HiCIAO dataset at a position angle offset by 
90 c with respect to the real disk. Figures 3d shows the re- 
sults of conservative LOCI applied to this modified dataset. 
The two disks appear virtually indistinguishable, demonstrat- 
ing that our simple model adequately explains the observed 
morphology. 


4.2. Surface brightness profiles 

We explore the viability of using the simple ADI image 
for the purpose of determining the II - band surface brightness 
profile of HR 4796 A’s debris disk. We extract a strip with a 
width of 1 1 pixels centered on the disk’s major axis from the 
intensity image, and collapse it along the minor axis. We re- 
peat this for the simulated disk image. The results are shown 
in a log-log plot in Figure 4. The 1 a error bounds are based 
on the radial noise profile in empty sky quadrants of the sim- 
ple ADI image before injecting the model disk. Simple ADI 
preserves the model surface brightness profile. The real disk’s 
observed profile is consistent with a single power law. 

5. DISCUSSION 

We have imaged and characterized the tapered outer bound- 
ary of HR 4796 A’s debris disk using ground-based angular 
differential imaging, and find that its iT-band surface bright- 
ness profile is well described by a simple power law with a 
slope around —9.5. Our observations are consistent with the 
expected result of a narrow planetesimal ring being ground 
up in a collisional cascade, yielding dust with a wide range of 
grain sizes. Radiation forces leave large grains in the ring and 
push smaller grains onto elliptical, or even hyperbolic trajec- 
tories. Simulations of such this process predict a smoothly 
tapered outer boundary, consistent with our findings (e.g., 
Krivov et al. 2006; Strubbe & Chiang 2006). Wahhaj et al. 
(2005) have demonstrated that the disk cannot be modeled 
with a single dust grain size, implying that dust populations 
of vatying grain size must be involved. The brightness slope 
is a result of changing dust densities and grain properties with 
distance. More detailed modeling is needed to relate the ob- 
served slope to the disk’s physical properties. 

Using our maximum regional merit technique, we have 
corroborated the evidence for an offset between the debris 
disk’s inner edge and the star. The dynamical influence 
of unseen planets orbiting within the disk cavity is most 
commonly invoked to explain such offsets (e.g., Kalas et al. 
2005; Thalmann et al. 2009; Buenzli et al. 2010, and refer- 
ences therein), and may be the cause of the proposed but un- 
confirmed intra-ring gap. Our confirmation of the ring offset 
also speaks against a dynamically cold source planetesimal 
population as proposed by Thebault & Wu (2008) as an alter- 
native to the planet hypothesis for explaining the morphology 
of the ring. The planets may well be detectable with the up- 
coming next-generation high-contrast imaging facilities. 

We thank David Lafreniere for generously providing us 
with the source code for his LOCI algorithm, and Jean- 
Charies Augereau for his GRaTer code. The authors acknowl- 
edge partial support from the Swiss National Science Foun- 
dation (SNSF), US National Science Foundation grants AST- 
1009203 and DGE-0646086, and a Japanese MEXT Grant-in- 
Aid for Specially Promoted Research (No. 22000005 ). 
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